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The locality of the fourth root of staggered fermion 
determinant in the interacting case. 
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The fourth root approximation in LQCD simulations with dynamical staggered fermions requires 
justification. We test its validity numerically in the interacting theory in a renormalization group 
framework. 
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1. Introduction 

Numerical simulations of dynamical fermions within the framework of the staggered formal- 
ism are both computationally cost effective and phenomenologically successful. Exploiting the 
advantageous properties of an improved staggered fermion formulation, various collaborations are 
performing high-precision lattice QCD calculations that are in excellent agreement with experi- 
mentally known measurements [jl|. However, this success is clouded by the long-standing prob- 
lematic issue of the validity of the fourth root approximation. The staggered fermion describes, in 
fact, four tastes on the lattice, so in order to study QCD with Nf = 2 or Nf = 2 + 1, a Boltzmann 
weighting is used that contains a fractional power of the fermion determinant. Because of taste- 
breaking at nonzero lattice spacing, taking the fractional power of the staggered determinant before 
restoring taste symmetry is conceptually nontrivial. The issue whether the fourth root prescription 
gives a lattice theory in the right universality class to reproduce QCD is yet unresolved and the 
phenomenological success of the staggered formulation has heated up the discussion. In the past 
years various numerical investigations [Q] have addressed this concern. Recently a direct approach 
was adopted towards establishing whether the universality class is the right one: if a single-taste 
local fermion action can be found, whose determinant is equal to the fourth root of the staggered 
fermion determinant up to cutoff effects 

lim(detD 4f ) 1/4 = detDdet// (1.1) 

(where D is a local Dirac operator and H is local and contains only cutoff effects) then the fourth 
root prescription can be consistently fitted into the framework of a local field theory. In the free 
theory numerical [^] and analytical studies [Q, |5|] showed that such a local operator exists. In 
particular, in ^ Shamir applied renormalization group blocking to the free staggered operator Do 
in the spin ® taste representation. After n blocking steps Q n , the fermionic degrees of freedom live 
on the coarse lattice with lattice spacing a c = 2"a and the determinant of the staggered operator 
decomposes as detD^a = 2~ n a c ) = detD,, detG^ 1 . More details can be found in [Q], where 
it was proved analytically that the blocked propagator D^ 1 = a -1 + Q n D^ factorizes, in the 
limit n — > °°, as D^ 1 = D r g ®/ and that G~ l = Do + a Q/ n Q n is a local operator that contains only 
cutoff effects. This completes the proof that in the continuum limit, i.e. a{= 2~"a c ) — ► 0, the 
decompositions of Eq. ( |L~l| ) holds in the free case. 

2. Renormalization group transformation in the interacting case 

In this section we briefly discuss the renormalization group (RG) program adopted to test 
the interacting case. The new complication is that a mapping from the one-component staggered 
fermion basis to the spin (8) taste basis is not unique. We define the interacting theory in the 
one-component formalism. The first fermion RG blocking transformation is used to perform the 
transition to the spin (g) taste representation: 

{Q {1)W ) aa {2x) = \ £ Wf 2 ^f 4 4 ) (aa) W r (2x,2 X + r) ¥ (2 X + r) 1 (2.1) 

^=0,1 
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where a and a are Dirac spin and taste indices respectively. We require the parallel transporters 
W r (2x, 2x + r) to be a sum over all the shortest paths 1 between the origin 2x and the other sixteen 
sites 2x + r of the hypercube. The sum over different paths is thought to reduce the breaking 
of hypercubic symmetry; however, the choice to parallel transport all the points within a given 
hypercube to the corner r = unavoidably introduces some asymmetries. The subsequent RG 
steps are a covariant generalization of the arithmetic mean on a 2 4 hypercube used in [||] : 



(0 (B) VO(2x)=2 



-4 



£ W r (2x,2x + r)\jf(2x + r) 

r u =0,l 



n > 1. 



(2.2) 



The shortest paths that compose the parallel transporters W r (2x,2x + r) in Eq. (2^2) are constructed 
from links of the blocked lattice (with lattice spacing 2" a). These blocked links are built following 
a program suggested in In detail, the links that live on the fine lattice (with lattice spacing 
2" _1 a) are twice APE smeared and projected back to SU(3). The blocked links are then built 
multiplying two of these fine smeared links in line. Once the blocking kernels are defined we 
project the nondiagonal piece of the blocked propagator D^ 1 onto the spin ® taste Clifford space: 



{QnD Q l Ql){x,y 



0)=£r s ®4iJf sr (je,y = 0). 

S,T 



(2.3) 



= 2 live on the coarsest lattice and Q n can be regarded as a 'big' blocking step 
i) ■ ■ ■ <2(i) that transforms the original lattice with lattice spacing a into the blocked 



Here x and y 

Qn = Q(n)Q(n 

lattice with lattice spacing a c = 2" a. The coefficients M$t {x,y = 0) of the projection in Eq. (2.3 ) are 
evaluated numerically and averaged over two different lattice ensembles. We fixed Lorentz gauge 
before doing any RG blocking. Parameters of the simulations are shown in Table |]. Applying a 
different number of RG steps, n, on each of these two ensembles, the resulting blocked lattices 
have the same coarse-lattice spacing of a c = 0.72 fm. This allowed us to see how the coefficients 
Mst scale with one additional RG blocking step. 



unblocked lattice 


16 3 x48 


40 3 x 96 


a 


0.18 fm 


0.09 fm 


m s a 


0.125 


0.05 


number of RGT's 


2 


3 


blocked lattice 


4 3 x 12 


5 3 x 12 


a c 


0.72 fm 


0.72 fm 


number of cfgs 


148 


56 



Table 1: Simulation parameters 



As a consistency check on the blocking procedure described above, we evaluated the mass 
of the (ss) meson. In fact, the spectrum on the blocked lattice is expected to be identical to the 
spectrum in the original theory. The meson mass can be evaluated directly from the coefficients 



M ST {x,y) ofEq. @: 



1 These are built using APE smeared links of the original lattice to be consistent with the RG program adopted for 
the subsequent RG steps. 

2 The fermion source is defined at the origin. 



114/3 



The fourth root approximation: interacting theory. 



F. Maresca 



cUt) = L < Tr(r 5 r r r r /rrT54,) \M sr (t,x)\ 2 >, (2.4) 

where the r$'s and the Ty r's are the 16 Dirac matrices in the normal and adjoint (t^ = y*) rep- 
resentation respectively. In Figure [l] we present the taste splittings for the meson masses evaluated 
on the two blocked lattices. At this level of statistics, the decrease in splitting from n = 2 to n = 3 
is consistent with the expected &{a 2 a) or 0(a 2 a 2 ), but the preliminary RG blocked splittings are 
systematically higher than values obtained in direct measurements at much lighter valence quark 
masses. Further study is needed. 
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Figure 1: Pseudoscalar taste splittings after 2 and 3 RG blockings in units of r\, the distance at which the 
static quark potential satisfies r 2 dV jdr = 1 . 



3. Preliminary results 

In this preliminary study we have analyzed for various displacements [jc — y| the dominant 



terms in Eq. ( |2.3| ) that are also present in the free theory: 



I (7m ®'M?W) + (75 <8>Wp)B$\x,y) + £ (75 ® x 5 x v )B' v {n \x,y) + (/®7) + 

(3.1) 

Comparing Figures |2] and |3[ it is remarkable to see that increasing by only one unit the number 
of RG blocking transformations the mass term (/ <g>I) ( □ ) becomes larger than the taste breaking 
term 75 ® TsT,- ( O ) for any direction i for displacement \x — y\ = 1. This is in agreement with what 
we expect from the free theory: in the limit n — > 00 the only terms in M^r that should survive are 
diagonal in taste space. 



Other terms in Eq. (|3JJ) are seen in our simulations. Their contribution is 10 — 1000 times 
smaller than the leading kinetic term in Eq. ( |3~l| ) and this make it difficult to identify the 
terms that are statistically relevant. We need to increase the number of configurations to have a 
better understanding of the statistics. Besides this issue, it is crucial to know how these terms 
scale with n to determine whether the blocked propagator becomes diagonal in taste space in the 
interacting theory as it does in the free theory. However an element that can play an important 
role in this scaling is the breaking of the hypercubic symmetry, caused by the introduction of the 
parallel transporters, as it has been briefly discussed in the previous section. In fact, Figures || 
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Displacement of 1 



16 x48 lattice 

n = 2; cfgs = 148 



- 4> <j, * 



~-FX +y +z +r 

(+) 



o Y| ® I 

o Y 5 ® TjT 

□ I®I 

a T 5 ® T 5 T k (k 5t i) 



=y =z =t~ 

(-) 



(1) 



Figure 2: M S7 - after n = 2 RG steps. 



Displacement of 1 



40" x 96 lattice 

n = 3; cfgs = 56 



O Yj ®I 

□ I®I 

o Y 5 ® iji, 

a Y 5 ® < k * ') 



* <> 4> * 

A ^ A A 

+x +y +z +T 

(+) 



(1) 



<t> <r> <t> 

4 A A A 

^x : y =2 T~ 
(-) 



Figure 3: M ST after n = 3 RG steps. 



and |3| suggest that the blocking procedure is introducing some hypercubic asymmetries, since the 
coefficients Msr{x,y) are not perfectly degenerate under reflections about the principal axes. In 
order to investigate further these terms, present only in the interacting theory, it is necessary either 
to quantify the hypercubic symmetry breaking or to redefine the RG blocking transformations in 
order to guarantee hypercubic invariance (one suggestion is presented at this conference 



4. Interacting theory vs free theory 

In ^ the scaling properties of the taste breaking terms By (p) were evaluated analytically. It 
was shown in momentum space that they scale like 2~ n , uniformly in p, so the blocked propagator 
becomes diagonal in taste space, when n — > oo. Postponing the issue of the breaking of hypercubic 
symmetry to a later study, we consider here a linear combination of the amplitudes |Msr(jc,;y)| 
under reflections and rotations about the principal axes. In Figures ^ and || we show how the ratios 
between the taste violating terms and the leading kinetic term A W scale with n. 
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Figure 4: Taste breaking term Z? (n) over the ki- Figure 5: Taste breaking term B /(n) over the ki- 
netic term. netic term. 

Our data show that they diminish when n is increased and that their magnitudes are also re- 
markably in agreement with the free theory represented by solid lines. In Figures ^ and ^ we show 
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the scaling properties of the taste breaking term divided by the mass term. The interacting theory 
agrees with the free theory for displacements \x — y\ = 1,2, however a discrepancy is seen at zero 
displacement where the interacting theory does not scale as expected. The zero displacement is not 
interesting when discussing locality thus this behavior, even if not well understood, does not spoil 
the good scaling DioDerties we are seeing in the interacting theory. 
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Figure 7: Taste breaking term divided by the mass 
term: free theory. 



Figure 6: Taste breaking term divided by the mass 
term: interacting theory 

The conclusion that we can draw from this preliminary study is that the blocked propagator in 
the interacting case is dominated by the same terms present in the free theory and that these terms 
scale as expected to make the decomposition of Eq. possible. The statistical relevance and the 
scaling properties of other terms seen in the interacting theory is still under investigation. 

We thank Yigal Shamir and Maarten Goltermann for helpful discussions during the confer- 
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